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I. Introduction 


Over the past few years, theoretical physicists have shown a great deal of interest in 
solving higher dimensional Schrodinger equation for various potentials[1-12]. Apart from 
regular well known potentials, also some complicated potential like ring shaped potential 
coupled with Coulomb potential, exponential potential, hyperbolic potential or their differ¬ 
ent combinations are also addressed by many authors [13-19]. 

Unfortunately, in this part of literature so far the solutions of higher dimensional Schrodinger 
equation with anharmonic potentials like combination of singular fractional power potential 
with Singular integer power potential V(r) = Y^=o Sextic potential V(r) = ar 6 + br A +cr 2 , 
Non-polynomial potential V(r) = r 2 + and others have not yet been paid a consid¬ 
erable amount of interest.Singular fractional power law potential has been recently found 
application in quark-antiquark interaction [20] and in models of exhibiting shape resonance 
behavior [21]. There are few number of studies basically on lower dimensional Schrodinger 
equation via various approximation method[22-25]. Motivated by the lack of these type of 
studies in higher dimensional Schrodinger equation, an attempt has been maid in this paper 
to solve the following square root power potential 


V(r) = ^ 

i =0 


cii 

T 5 

T 2 


( 1 ) 


where cq(i = 0,1, 2, 3,4) represent the potential parameters. The potential under consid¬ 
eration is not purely fractional power potential rather than it is a combination of singular 
fractional power potential with singular integer power potential. It is clear clue to the 
mathematical difficulty only the running index is taken up to i = 4. At present there are 
number of suitable techniques available [26-29] to solve different types of differential equation 
that generally appears in physical systems like Schrodinger equation in quantum mechan¬ 
ics. However, in this paper traditional wave function ansatz concept [30-31] with a little 
series solution technique [32] has been employed to solve the potential for multidimensional 
Schrodinger equation. 

The paper is organized as follows: Next section is for the bound state spectrum of the po¬ 
tential. In Section III special cases are studied. Section IV is for the conclusion of the work 
and finally important references are given in section V. 
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II. Bound State Spectrum of the Potential 


The D-dimensional time-independent Schrodinger equation for a particle of mass M with 
arbitrary angular quantum number £ is given by [33] (in units of h — 1) 

( 2 ) 


V D + 2 M yE ni - V (r) J J iftnim (j~ i CLd) = 0 , 

where E n £ and V(r) denote the energy eigenvalues and potential. Qp within the argument of 

n-th state eigenfunctions ^ denotes angular variables 61,62, 03, . ,&d- 2,P- The Laplacian 

operator in hyperspherical coordinates is written as 


Vr> = 


1 d D _ x d 


D r D -1 Q r 


dr' 


( 3 ) 


where 


A 2 - 


ST _I_ x ( 1 d cin k-i e d 

^ sin 2 6k+isin 2 6).+2 . sin 2 9 D _ 2 sin 2 p \sin k ~ l 9 k d9 k d9 k 

1 d . D _ 2 d 

+ . D _2 7T sm ^ o - 

sin u A p dp dp 


( 4 ) 


A|,_ 1 is known as hyperangular momentum operator. 

Expressing the solution as V-Vm( r - ^d) = -Rn7(' r )b7"(^Li) and applying the separation vari¬ 
able method in Eq.(2) we can have two separated equations 


A I^YP^d) = £(£ + D- 2)Y e m (Q D ), 


( 5 ) 


and 


d 2 D-ld £(£+D- 2) 2 Mr Tr/N 

or- + ~l~Tr r! + ~»\ Ent ~ v(r) 


Rnt{r) = 0 , 


( 6 ) 


where Y^ m ( fi#) in the Eq.(5) are the hyperspherical harmonics. £ is called the orbital angular 
momentum quantum number, where £ = 0,1, 2, 3. 

Eq.(6) is called the hyperradial or in short “radial” equation containing the radial part of 
the solution R n i{r). The separation constant £{l + D — 2) is only applicable for D > 1. 
Now inserting the Eq.(l) into Eq.(6) and changing the variable as = x we have 


d 2 2D — 3 d 
dx 2 x dx 


N d 

X 2 


+ 8 M | ( E n f — a 0 )x 2 — aix — a 2 —— | 


R(x) = 0 , (7) 
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where = 4£(£ + D — 2) + 8Ma 4 and R n e(r) -H- i?(x). 

In order to solve Eq.(7) we first find the asymptotic behavior of R(x). As x —» oo we have 

- ~I~P~ + 8M j(£ n ^ - a 0 )a; 2 - aix|i?(x) = 0 . (8) 

The solution must be bounded at infinity i.e f?(x —» oo) —> 0. So let us try the solution of 
Eq.(8) as 

f?(x) = [const ] e ~( Ax2 + Bx ) ) (9) 


where A, B are constants. Taking this form of solution in Eq.(8) and neglecting the term 
B 2 — 2A, as it is fare for large limit of x, we have by comparing the coefficient of x 2 and x 


A — \/—2M(E n £ — ao), 

2 Mai 

15 = 


(10a) 

(10b) 


a/— 2M(E n £ — ao) 

The solution must be bounded in origin also. The well-behaved complete solution of R(x 
now can be taken as 


R(x) = [const]e- {Ax2+Bx) f(x ), 


( 11 ) 


were f(x —> 0) —» 0. 

Using Eq.(ll), Eq.(7) leads to the following homogeneous linear second order differential 
equation 

d 2 f(x ) ^df(x) 

-^ + hix) te +9{x)f{x) = °' (12) 

where 


h(x) 

g(x) 


2D-3 

x 
>2 


(. B 2 - 2 A- 8Ma 2 ) 


(4 Ax + 2 B), 

8Ma 3 + (2D - 3) (2 Ax + B) 


x 


N d 

iV ( 0,4 


X A 


(13a) 

(13b) 


The complete solution of this differential equation can be found if one integral in its com¬ 
plementary function be known. Let us try an inspection such that f(x) = x k is a solution 
of this equation. So we have 


k(k — 1) + H(x)k + G(x) = 0 , 


(14) 
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where H{x) = xh{x) and G(x) = x 2 g(x). 

Now h(x) and g(x ) are not analytic near the origin, but H(x) and G(x) are both analytic 
at x = O.This provides k(k — 1) + H( 0)k + G(0 ) = 0 and hence 

k+ = —(D - 2) + y/(D + 27 - 2) 2 + 8Mfl 4 , (15a) 

k. _ = -(D - 2) - a/ (D + 27 - 2) 2 + 8 Ma A . (15b) 


The acceptable solution must be k + (> 0) given by Eq.(15a) because it preserves the rule 
f(x -> 0) 0. 

Now we are in a situation to write the form of complete solution for f(x) as a series like 

OO 

f(x) = x k ^ c p x p . (16) 

p=0 

Taking this into Eq.(12) we have the following identity 

OO OO OO 

£ a p c p x k+p+2 - ^2 /3pC p x fc+p+1 + ^2 lpC P x k+p = 0 , (17) 

p =0 p =0 p =0 

where 


a p = B 2 — 8 M a 2 — 471(77 + k + p — 1), (18a) 

/5 P = 8Ma 3 + B(2D + 2/c + 2p — 3), (18b) 

7 P = (A: + p)(fc + p + 2D-4)-iV^ 4 . (18c) 


Equating the coefficient of lowest power of x, which corresponds to p = 0 in the left hand 
side of the identity, we have the indicial equation 

<*[*(*-1) + fc(2D - 3) - = 0. (19) 

We must choose Co ^ 0 as k = k + makes the bracket term of Eq.(19) a vanishing. 

Similarly equating the coefficient of next higher power of x. i.e x k+1 from the identity we get 

_ 8Ma 3 + B(2k + 2D — 3) 

Cl “ (*+l)(fc + 2D-3)-< 1 Co ' 

After deciding the coefficients Co and ci, finally equating the coefficient of x p+k+2 from both 
sides of the given identity i.e Eq.(17) we have 


OipCp (5pj r \c p j r \ T 7p-|_2Cp_)_2 0 . 


( 21 ) 
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Now for physical cases the series must be a polynomial or in short it must terminate at 
p = n. That means we must have c n ^ 0 and c n+ 1 , c n+2 , c n+3 ,... = 0 

From the above argument it is easy to find a n = 0 which finally provides the energy eigen¬ 
value equation in association with Eq.(lOb) as 

4 A 3 (D + k + n - 1) + 8 A 2 Ma 2 - 4 M 2 a\ = 0 . (22) 


The corresponding radial eigenfunctions are 


Rrd(r) 


N n (.e 


r^/-2M(E nl -a 0 )+r 1/2 


2Ma| 

\Z-2M(E ne -a 0 ) 


k± 

f 2 


n 



3=0 


(23) 


where we have used Eq.(10a, 10b),Eq.(ll) and Eq.(16) with the exchange of variable 
x = r 1 / 2 . N n £ is the normalization constant which can be evaluated from the condition 
Jo°[Rne( r )] 2 r D ~ 1 dr = 1 . 


III. Discussion of Special Cases 


1. Singular One Fractional Power Potential: V(r) = + ^§72 

The potential is obtained from Eq.(l) by setting a 0 = a 2 = a 4 = 0 
Under this condition from Eq.(15a) we have k + = 2£, and from 
Eq.(lOa) A = y/—2ME n e. 

Now Eq.(22) immediately gives the energy eigenvalues 


1 ( M 2 a\ V 
nl 2M\2k' + n+l) ’ 

where k' = (y + £ —■ 1 ). 

The corresponding radial eigenfunctions are given by 


Rne(f ) Al n (.£- 


rv^MUi+r 1 / 2 


y/-2ME n 


3 =0 


Cj r 


3 /2 


The above results are in excellent agreement with the ref. [33]. 


2. Mie-type Potential: V(r) — ao + y- + f? 

The potential is obtained from Eq.(l) by setting oq = a 3 = 0. 
Eq.(lOb) gives B = 0 hence Eq.(18b) provides f3 p (or /?)= 0. 


(24) 


(25) 
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Under this condition Eq.(20) gives C\ — 0 and the recurrence relation comes out from 
Eq.(21) as 


OlpCp T r )/p+2Cp+2 0 . 


(26) 


The above relation immediately gives 


{ =0, if p is odd, 
7 ^ 0, if p is even. 


(27) 


Taking p = 2 n r where n r = 0,1, 2, 3.we again argue that the series must terminate 

at p = 2 n r (i.e C 2 n r ^ 0) to achieve a polynomial solution to explain the physical cases. 
So under this situation in Eq.(26) we have to set c p+2 = 0 for p > 2 n r . 

This gives a 2 n r = 0 which is nothing but the same as Eq.(22) except p=n — 2 n r and 
potential parameter esq = 0. 

Hence the energy eigenvalues are 


En T £ Oq 


M 


0-2 


n r + l\l + v / (T) + 2£-2) 2 + 8Ma 4 


(28) 


This is what obtained in ref. [12] and ref. [31]. 

The radial eigenfunctions for this case emerges from Eq.(23) as 




N n «e 


-2M(E nr (-a 0 ) 


H 
y* 2 


2 fi}’ 

c 1 r J ' /2 . 


j =0 
j^odd 


(29) 


3. Coulomb Potential: V[r) = 

This can be obtained from Eq.(l) by setting a 0 = (i\ = o 3 = a 4 = 0. Thus the results 
of Mie-type potential as well as the explanations of that section are equally applicable 
and can be used here with an extra condition a 0 = a 4 = 0. So the energy eigenvalues 
are 


En r i 


a 2 


M 


2 \ y + £ + ti. 


_ 1 

r 2 


(30) 


This result is similar to the ref. [10]. 

Eq.(15a) provides k + = 2£, hence we have the radial eigenfunctions from Eq.(29) as 




—r 

Nn r £& 


yj 2 ME nr £ 


E 


Cj r 


,3 /2 


j=0 

j^odd 


(31) 
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IV. Conclusions 


In this paper D -dimensional Schrodinger equation is addressed for square root power law 
potential using wave function ansatz method. It has been shown that traditional mathe¬ 
matical method like wave function ansatz and series solution always provide a wide space to 
deal with complicated potential problems. The result that has been obtained in this work is 
satisfactory in the sense that some special cases like singular one fractional power potential, 
Mie-type potential, Coulomb potential are derived easily and they are in excellent agreement 
with the previous works. 

The exact form of the energy eigenvalues given by Eq.(22) is complicated. The only way 
to find the energy eigenvalues is the numerical computation of the equation for a particular 
state. However the eigenfunctions are in well known form. May be after normalization they 
can be expressed in terms of Laguerre polynomials or confluent hypergeometric functions. 
We will look forward to the further study of this potential in future at least for few more 
terms {i > 4). Any other method like asymptotic iteration method (AIM), Nikiforov-Uvarov 
method (NU) will be helpful to this. 

Furthermore, following the same way of this paper, one can study the other class of fractional 
singular power potentials like ^ and ^ . 
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